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1.  INTRODUCTION 


We  consider  model  fitting  problems  in  which  the  data  are  fuzzy  vectors, 
and  the  mathematical  model  of  the  observed  event  is  a  system  of  equations 
involving  observables  and  mode!  parameters.  Such  problems  naturally  arise 
when  data  are  scarce  and  information  is  not  available  about  distributions  of  data 
inaccuracies.  Then  it  may  only  be  possible  to  obtain  estimates  of  data  scatter 
from  past  experience  in  the  form  of  estimated  data  membership  functions. 

The  fitting  of  the  model  consists  of  finding  an  optimal  set  of  parameters 
and  a  corresponding  set  of  corrections  to  the  observations  such  that  the  model 
equations  are  satisfied.  One  obtains  different  solutions,  depending  on  the 
definition  of  optimality.  In  this  article,  the  optimally  is  defined  in  terms  of  a 
least  squares  objective  function.  Particularly,  we  shall  maximizing  the  sum  of 
squares  of  the  deviations  of  the  membership  values  of  one.  One  obtains  rela¬ 
tively  simple  algorithms  for  this  problem  if  the  membership  functions  of  the 
data  vectors  belong  to  a  particular  class  of  functions.  This  class  of  conical  func¬ 
tions  is  defined  and  discussed  in  Section  3,  after  a  short  exposition  of  general 
fuzzy  model  fitting  problems  in  Section  2.  The  model  fitting  to  fuzzy  vectors 
with  conical  membership  functions  is  described  in  Section  4,  and  Section  5 
presents  an  application  example. 

An  obvious  extension  of  the  problem  to  fuzzy  models  is  mentioned  in  Sec¬ 
tion  2,  and  will  be  treated  in  a  forthcoming  paper. 


2.  FUZZIFICATION  OF  MODEL  FITTING  PROBLEMS 


We  define  in  this  section  a  general  least  squares  model  fitting  problem  and 
outline  several  possibilities  how  the  problem  can  be  fuzzified. 

Let  XfEn  be  the  vector  of  all  observations  of  an  event,  and  let  the 
mathematical  model  of  the  event  be  the  equation 

FIX, t)  =  0  (2.1) 


1 


where  F  is  a  8-dimensional  function,  and  t  is  a  p-dimensional  model  parame¬ 
ter.  We  assume  that  F  is  twice  differentiable  with  respect  to  all  its  arguments, 
and 


n  >  s  >  p  >  0  (2,2) 

The  goal  of  a  model  fitting  operation  is  to  find  such  a  parameter  vector  t  that 
the  model  equation  (2.1)  is  satisfied.  Generally,  this  goal  cannot  be  achieved, 
because  the  number  s  of  equations  in  (2.1)  usually  is  larger  than  the  number  p 
of  free  parameters.  Then  a  solution  of  the  problem  only  can  be  obtained  by 
weakening  the  requirement  that  Eq.  (2.1)  be  exactly  satisfied.  The  replacement 
of  Eq.  (2.1)  can  be  done  in  different  ways  and  one  obtains  corresponding  dif¬ 
ferent  solutions. 

In  usual  least  squares  model  fitting  one  assumes  that  the  observations  X 
may  contain  errors  and,  therefore,  the  model  equations  should  not  be  satisfied 
at  A",  but  in  the  vicinity  of  X.  That  is,  one  corrects  the  observations  A’  by  a 
term  cy  and  replaces  the  model  equations  (2.1)  by  the  constraint  equations 

F{X+cx,t  )  =  0  (2,3) 


One  would  like  to  have  only  small  corrections  and  this  is  achieved  by  seeking 
such  a  solution  cx  and  t  of  Eq.  (2.3)  which  minimizes  the  square  of  a  norm  of 
fy.  The  problem  may  be  formulated  as  follows: 

minimize  W  =  ||  J| 

subject  to  F(X+  Cx  ,  t  )  =0 

Hence,  a  least  squares  model  fitting  problem  is  defined  by  the  following  set 

{  X  •  II  CX  II  i  F  i  W  -  ||  cx  ||*  }  .  12-5) 

containing  a  data  vector,  a  definition  of  a  norm,  a  definition  of  a  model  func¬ 
tion,  and  a  definition  of  an  objective  function  as  the  square  of  the  norm  of 
corrections. 

The  problem  (2.5)  can  be  fuzzified  as  follows.  First,  the  data  A'  may  be 
fuzzy.  In  this  case  one  has  in  addition  to  the  observed  X  also  a  definition  of  its 
membership  function  p  x.  Then  we  replace  the  definition  (2.5)  of  the  least 
squares  model  fitting  problem  by 


(2,4) 
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{x  ;  /ix(cx)i  Fi  W=(1-MX(CX))2]-  ,  (2.6) 

or 

(2«7) 

subject  to  F  (  X  +  Cx  ,  t  )  s  0  .  J 

This  type  of  model  fitting  problems  with  fuzzy  data  is  treated  in  Section  4. 

Another  possible  fuzzification  of  the  problem  (2.5)  is  obtained  by  intro¬ 
ducing  a  fuzzy  model  function.  This  can  be  done,  for  instance,  by  defining  a 
set  of  functions  with  an  associated  membership  function  / if  The  data  X 

then  are  assumed  to  be  crisp,  that  is,  one  assumes  <^=0.  The  corresponding 
least  squares  model  fitting  problem  is  defined  by  the  set 

{x  ;  ||cx||  =  0  ;  F;  MF  ;  vv=  (l-/ip)2}  ,  (2.8) 


minimi ze 


W  =  (  1-M)(  ( Cx) ) 


or 


minimize 
subject  to 


W  =  (l-/xF)2 
F  (  X ,  t  )  =  0  . 


(2.9) 


Problems  of  this  type  are  akin  to  a  problem  formulated  by  Tanaka  et  al.  We 
shall  treat  the  problem  (2.8)  in  a  forthcoming  paper. 

Another  kind  of  fuzzification  is  a  combination  of  the  former  two,  that  is, 
one  assumes  that  the  data  are  fuzzy  and  use  a  fuzzy  function  to  fit  the  data.  It 
is  probably  the  most  useful  case  for  applications. 

Finally,  one  may  fuzzify  the  objective  function  VP,  for  instance,  by  allow¬ 
ing  to  minimize  other  powers  than  the  square  of  the  norm  or  membership 
function.  Because  we  restrict  the  present  investigation  to  least  squares  minimi¬ 
zation,  we  shall  not  further  discuss  this  typf  of  fuzzification. 


3.  FUZZY  VECTOR  SPACES  AND  SPREAD  PROPAGATION  FORMULAS 


In  many  model  fitting  problems  the  data  consist  of  a  set  of  vectors  rather 
than  a  set  of  scalars.  Such  data  can  be  treated  in  the  context  of  fuzzy  sets  by 
defining  for  each  vector  a  membership  function  that  depends  on  all 
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components  of  the  vector.  By  chosing  a  proper  class  of  membership  functions 
one  can  simplify  the  numerical  handling  of  the  corresponding  model  fitting 
problem.  We  shall  introduce  and  discuss  in  this  section  such  a  class  of  func¬ 
tions,  namely  a  particular  extension  of  the  triangular  membership  function  of  a 
fuzzy  scalar  to  n  dimensions. 

Let  a  be  a  fuzzy  scalar.  We  associate  with  it  the  triangular  membership 
function  pt(x),  defined  by 


Ma  ( x  ) 


max 


(3.1) 


We  call  a  the  apex  and  st  the  spread  of  a  and  notice  that  pa(x)  is  completely 
determined  by  the  two  quantities  a  and  «a.  The  spread  may  be  considered  as  a 
measure  for  the  accuracy  of  a.  Then  a  reasonable  dimensionless  measure  for 
the  difference  between  a  crisp  point  x  and  the  fuzzy  observation  a  is  the  norm 


x-cf  =  x-a/s. 

a  1 


(3.2) 


The  relation  between  this  norm  and  the  membership  function  of  a  is 


Mn<X)  =  1-  min  jl  ,  \\x-7i\\Q  j  .  (3.3) 


Now  we  extend  this  simple  function  to  n  dimensions  and  demand  that  the 
extension,  that  is,  the  membership  function  ftA(  X)  of  a  fuzzy  vector  .4,  has 
the  following  properties: 

(A)  Each  component  of  .4  has  a  triangular  membership  function,  and 

( B)  the  function  ( l-/i  4 ( A') )2  is  differentiable  in  the  domain  where 
/^(A')>0. 

The  latter  property  is  important  for  the  simplification  of  the  numerical  treat¬ 
ment  of  model  fitting  problems. 

A  function  with  the  properties  (A)  and  (B)  is  the  conical  membership 
function 

Ma(X)=  max  jo,  1-  [(X-A)T  P'Mx-A)]^2  ,  (3.4) 

where  PA  is  diagonal  matrix  with  diagonal  elements  which  are  the  spread 
squares  sj;  of  the  components  d,  of  the  vector  .4  .  Figure  1  shows  pA(X)  in 
the  case  n  =  2. 
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In  analog)'  to  Eq.  (3.2)  we  also  define  a  norm 

llX_^llA=  [(X-A^P^tX-A)]172  (3.5) 

and  obtain  an  expression  of  the  membership  function  in  terms  of  the  norm: 

^A(X)  =  1  -min-|l ,  ||  X-A  jj^  .  (3.6) 

The  boundary  of  the  support  of  the  membership  function  is  given  by  the 
equation 

(X-A)T  P’Mx-A)  =  1  .  <3<7) 

This  boundary  is  a  hyperellipsoid  in  the  n-dimensional  X-space.  The  principal 
axes  of  the  ellipsoid  are  parallel  to  the  coordinate  axes,  and  their  lengths  are 
equal  to  2sai-.  The  matrix  PA  which  defines  the  membership  function,  we  call 
the  panderanee  matrix  of  the  fuzzy  vector  A. 

It  is  practical  to  slightly  generalize  the  given  definition  of  a  conical 
membership  function  by  permitting  the  support  ellipsoid  to  have  principal  axes 
that  are  not  parallel  to  the  coordinate  axes.  Such  a  situation  arises,  e.g.,  if  one 
choses  a  different  coordinate  system.  We  determine  the  proper  generalization 
of  the  definition  by  considering  a  linear  coordinate  transformation 

2  =  DX  +  Z0  .  (3.8) 

The  apex  .4  of  the  fuzzy  vector  .4  with  the  panderance  matrix  PA  is  mapped  by 
this  transformation  into  a  vector 

B  =  D A  +  Z0  ,  (3.9) 

and  the  panderance  matrix  Pp  of  D  is  obtained  by  the  following  spread  propa¬ 
gation  formula: 

PB  =  DPADT  *  (3. 10) 

We  first  verify  this  formula  for  the  ca.n-  where  (3.8)  is  a  regular  coordinate 
transformation,  that  is,  D  is  a  non-singular  (  n  Xn) -matrix.  Then  Pg,  computed 
by  Eq  (3.10),  is  positive  definite  and  the  membership  function 

Me(Z)  =  1-minjl  ,[<Z-B)T  p-1(Z-B)]l/2J  (3.11) 

has  the  property 


^tg(Z  (X) )  -  /aa  (X  )  . 


(3.12) 
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because 


(3.13) 


(Z-B)TPgl(Z-B)  =  (X-A)TDT(DT)-1Pa-1D'1D(X-A)  =  (X-A)TPAl(X-A) . 

In  the  Z-cocrdinate  system,  each  component  again  has  a  triangular  membership 
function,  the  projection  of  the  cone  defined  by  the  panderance  matrix  PB.  The 
spread  sbi  of  each  component  i>,  of  B  equals  the  square  root  of  the  i-th  diago¬ 
nal  element  of  PB.  The  support  ellipsoid’s  principal  axes  are  in  general  not 
parallel  to  the  axes  of  the  Z-coordinate  system. 

One  sees  from  the  formula  (3.10)  that  the  panderance  matrix  of  a  conical 
membership  function  merely  has  to  be  positive  definite.  However,  as  we  shall 
see  later,  degenerated  support  ellipsoids  do  occur.  Therefore  we  shall  assume 
that  the  panderance  matrix  only  is  positive  semidefinite. 

The  elements  of  a  panderance  matrix  PA  have,  in  general,  different  physi¬ 
cal  dimensions,  corresponding  to  the  dimensions  of  the  components  of  A . 
They  also  may  have  different  orders  of  magnitude,  depending  on  the  units 
which  are  used  to  express  A.  Therefore,  a  more  lucid  representation  of  PA  is 
in  terms  of  the  spreads  sa-  of  the  components  a-  of  A  and  a  concordance  matrix 
C'A,  which  is  related  to  PA  by  the  formula 


or 


PA  s  SA  CA  SA  » 


-1 


PASA 


(3  .14) 


where  SA  is  a  diagonal  matrix  with  spreads  sai  in  the  diagonal.  Elementwise 
the  relation  (3.14)  between  panderances,  concordances  and  spreads  is 

A  i  j  ^Aij  ^oi^Qj’  '  »  1  -  1  »  .  . .  ,  n  .  (3.15) 

The  concordance  matrix  has  dimensionless  elements,  less  or  equal  one  in  abso¬ 
lute  value,  with  ones  in  the  diagonal. 

So  far  we  only  have  verified  the  spread  propagation  formula  (3.10)  for  reg¬ 
ular  linear  coordinate  transformations.  It  is,  however,  also  valid  for  any  other 
linear  function  (3.8)  in  the  sense  that  the  membership  function  which  is  pro¬ 
duced  by  Eq.  (3.10)  for  each  component  of  Z  satisfies  the  extension  principle. 
Now  we  shall  prove  this  property  of  the  formula. 

Let  one  component  of  the  linear  function  (3.8)  be 

Z  =  dTX+z0  .  (3.16) 
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The  fuzzy  vector  .4  is  mapped  into  a  fuzzy  number  b  with  the  apex 

b  =  dTA  +  z0  .  (3.17) 

The  panderance  (the  spread  square)  of  b  is  according  to  (3.10) 

sb2=dTPAd  (3.18) 

and  the  membership  function  of  fc  is 

Mb(z)  *  1- min(l  ,  |  z-b|/»b)  .  (3.19) 

According  to  the  extension  principle,  the  membership  function  of  should  be 

Mb(z)  *  max  ^(X)  •  (3.20) 

X :  dTX  +  z0  =  z 

We  shall  show  that  /ij,(r),  computed  from  (3.20),  equals  the  computed  with 
the  formula  (3.10).  To  that  end,  we  reformulate  the  right  hand  side  of  Eq. 
(3.20)  as  the  following  constrained  minimization  problem: 

minimize  ) j  X -A  ||  A  =  (X-A)"^  ^(X-A) 

subject  to  d^(X -  A)  =  z-b  . 


We  use  a  Lagrange  multiplier  k  and  obtain  the  following  normal  equations  of 
the  problem  (3.21): 

Pa_1(X-A)  - kd  =  0  , 
dT(X-A)=  z-b  . 

We  solve  the  first  equation  (3.22)  for  k,  obtaining 

k  =  dT(X-A)/(dTPAd) 

or 

k  =  (z-b)/(dT^d)  .  (3.23) 

Substituting  this  into  the  first  equation  (3.22)  and  multiplying  from  the  left  by 
( A  -  A  )  T  one  obtains 


o  r 


llx_sliA  min=  l-»K  • 


(3.24) 

(3.25) 
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Therefore 


mox'V*  )  •  1  -  minjl  ,||x-A||A  mJ  -  l-min{l,|z-b|  /.„} .  (3.26) 

that  is,  the  solution  of  (3.20)  is  (3.19). 

The  proof  can  be  modified  to  cover  also  singular  cases  with  semidefinite 
panderance  matrices  PA  by  restricting  the  minimization  (3.21)  to  the  subspace 
where  PA  is  definite.  Hence,  we  have  show’n  that  the  spread  propagation  for¬ 
mula  (3.10)  componentwise  satisfies  the  extension  principle  for  any  linear  func¬ 
tion  of  X.  If  the  function  Z(X)  is  not  linear,  then  one  may  obtain  an  approxi¬ 
mate  panderance  matrix  PB  for  B=Z(A)  by  linearizing  the  function.  Then 
the  propagation  formula  is 

p  .  az  D  /az\T 

p8  •  ax  PAVax)  •  <3.27) 


where  dZ  jd X  is  the  Jacobian  matrix  of  the  function  Z(  X).  The  quality  of  the 
approximation  depends  on  the  size  of  the  second  order  derivatives  of  Z ( A") , 
relative  to  the  panderances  and  to  the  first  order  derivatives. 

So  far  we  only  have  considered  single  fuzzy  vectors  and  defined  a  norm 
for  the  distance  between  a  fuzzy  and  a  crisp  vector.  The  membership  function 
was  expressed  by  Eq.  (3.6)  in  terms  of  this  distance.  We  shall  also  need  a 
measure  for  the  distance  between  two  fuzzy  vectors.  Now  we  give  the  neces¬ 
sary  definitions  and  take  for  simplicity  a  general  approach. 

Let  Hn  be  a  vector  space  with  fuzzy  n-component  vectors  as  elements. 
We  assume  that  each  A(Hn  is  assigned  a  metric  which  defines  a  distance 
A(X,.4)  in  /?„,  that  is,  between  any  X(Rn  and  AcHn.  We  also  assume  that  the 
membership  function  ^a[X)  of  .4  is  related  to  h(X,A)  by  the  formula 

A*A(X)  =  max  |o,  l-h(X.A)}'  =  1-  min|l,  h(X,X)J  .  (3.28) 

The  distances  h(X,A)  are  not  useful  for  the  measurement  of  distances  in  Hn 
because  in  general  h[A,B)^  h(B,A).  We  therefore,  define  a  discord  between 
two  elements  of  //„  bv  the  formula 


D(  A,  B)  =  min  max 
X«FL 


|h(X,A),  h(X,B)j  (3.29) 


D  is  symmetric  but  it  is  not  a  distance  in  Hn,  that  is,  it  has  the  properties 


ft 


and 


D{ A,  B)  >  0  if  A#B 
D(X,§)  *  0  if  A  =  B 


(3.30) 


but  it  does  not  satisfy  the  triangle  inequality. 

As  a  complement  to  the  discord  D(A ,  B)  we  also  define  a  grade  of  colloca¬ 
tion  of  two  fuzzy  vectors  A  and  B  by 

y(A  ,  B  )  *  max  min  jfiA(X), /ig(X)j>  .  (3.31) 

x*  Rn 

Because  we  have  assumed  the  relation  (3.28)  between  ft A  and  h(A',A),  we  also 
can  express  the  grade  of  collocation  in  terms  of  the  discord: 


y  (A,B)  =  max  jo,  1-D(X,§)|*  1-  min|l ,  0(A,B)|  .  (3.32) 


In  the  special  case  of  vectors  with  conical  membership  functions  we  define 
h  ( A\  A  )  =  |  |A-  A  1 1.  In  th  is  case,  the  distance  h(X,A)  is  inversely  proportional 
to  the  spreads  of  the  components  of  A.  If  one  multiplies  the  spreads  of  B  by  a 
factor  d>  and  lets  <p  approach  zero,  then  one  obtains  for  the  discord  the  limit 

lim  D(B,  A)  =  h(B,  A)  .  (3.33) 

<£-•“0 

The  limit  of  a  membership  function  with  vanishing  support  may  be  considered 
as  the  membership  function  of  a  crisp  vector.  In  this  sense,  the  discord 
D[B,A)  approaches  the  di  tance  h(B,A)  between  a  crisp  vector  and  A,  as  B 
becomes  a  crisp  B.  Hence  the  discord  is  a  convenient  generalization  of  the  dis¬ 
tance  in  Rn  with  a  proper  limit  property. 
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4  LEAST  SQUARES  MODEL  FITTING  TO  FUZZY  DATA 


A  general  formulation  of  a  least  squares  model  fitting  problem  with  fuzzy 
data  was  given  in  Section  2  by  Eq.  (2.7).  Now  we  shall  give  &  slightly  different 
formulation  which  is  more  adequate  for  the  present  discussions. 

Let  the  observations  be  fuzzy  vectors  Xit  i=l, each  having  the 
dimension  n,,  and  a  membership  function  fijg,  represented  by  a  panderance 
matrix  Px.  Let  the  model  of  the  event  by  defined  by  the  equations 

Fj  (Xj,  t)  =  0,  i=l,...,S,  (4.1  ) 

where  t  is  a  p-dimensional  model  parameter  and  F ’•  are  r, -dimensional  func¬ 
tions.  We  assume  that  the  dimensions  n,  r,,  and  p  satisfy  the  conditions 

I" i  >  p  >  o  .  (4.2) 

1  1 

Compared  to  the  formulation  in  Section  2  we  now  have  assumed  that  the 
model  equations  have  subsets  which  depend  on  distinct  subsets  of  the  observa¬ 
tions.  (The  present  formulation  includes  the  formulation  of  Section  2  as  the 
special  case  5=1)  c\~h  a  partitioning  is  not  always  possible,  but  it  offers  sub¬ 
stantial  algorithmic  advantages.  I3!  Also,  many  model  fitting  problems  are 
naturally  formulated  by  partitioned  equations  (4.1),  e.g.,  in  standard  least 
squares  problems  where  the  Ft  are  scalar  functions  and  all  A '  are  distinct. 

A  least  squares  model  fitting  problem  now  may  be  formulated  as  follows: 


minimize  W  =  £  [l-/xx-  (Xj  -•*  ep]  , 
subject  to  Fj  ( X;  +  Cj  ,  t )  =  0  ,  i  *  1 ,  .  .  .  ,  $  , 


(4.3) 


where 


MXj(Xj  +  cj  )  =  1-m 


in|lI(cjTP|r1c.)1/2J  .  (4.4) 


10 


We  recall  the  definition  (3.5)  of  a  norm 


(4.5) 


and  observe  that  the  corrected  observation  A^-f-c,-  only  is  meaningful  if  the 
correction  r,  satisfies  the  inequality  ||c,|j^7<  1-  If  lkill\«>  1  then  the 
membership  function  ^.\7(-V,  +  c|)  vanishes  and  in  general  the  function 
Fj(A,+e(  t)  might  not  even  be  defined.  Therefore,  one  should  supplement  the 
constraints  of  the  minimization  problem  with  the  conditions  |  |f,  |  |.\7  <  l-  1  with 
a  satisfying  0<-j  <1. 

We  also  express  the  objective  function  in  terms  of  the  ||e,  ||.\72  using  the 
relation  (4  -4).  The  advantage  of  this  reformulation  is  that  the  new  objective 
function  is  differentiable  with  respect  to  the  c,  for  all  finite  ct.  As  a  result  one 
can  use  Lagrange  multiplier  techniques  for  the  numerical  solution  of  the  prob¬ 
lem.  The  restriction  to  conical  membership  functions  was  made  in  order  to 
have  such  a  differentiable  objective  function.  Using  the  new  objective  function 
and  the  above  discussed  additional  constraints  one  obtains  the  following  con¬ 
strained  minimization  problem: 


minini ze 

W  =  £  Cj  Pj  Cj  , 

1  [ 

(4.6) 

subject  to 

Fj(  Xj+cj  ,t)  =0,  i*i,  ...,s  J 

and 

||  ci || xj  <  •  i *  1 .  - »  • 

(4.7) 

One  observes  that  the  formal  difference  between  an  ordinary  least  squares 
problem  and  a  model  fitting  with  fuzzy  data  are  the  constraints  (4.7).  Concep¬ 
tually,  these  constraints  are  compatibility  conditions  between  data  and  model, 
and  they  identify  outliers.  Geometrically,  they  require  that  for  all  i  =  l,...,s  the 
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ill  m 

support  of  the  function  ftx{X)-  7  intersects  with  the  fitted  surface  /,( X,t)—0. 
However,  the  results  of  the  fitting  is  not  a  crisp  function,  but  a  fuzzy  set  of 
functions,  because  it  has  been  computed  from  fuzzy  data.  Therefore,  the  fitted 
surface  also  is  fuzzy.  In  the  A'-space  it  has  a  membership  function  ftFi(X), 
which  equals  one  on  the  surface  F,  =  0,  and  has  a  finite  support  around  that 
surface.  The  condition  (4.7)  means  that  the  support  of  the  functions 
/z AT( A')  —  *  and  /ifl(A')-l  are  required  to  intersect.  It  is  more  reasonable  to 
only  require  that  the  supports  of  Vx;(X)~  1  and  q  intersect.  Now  we 

shall  formulate  a  condition  based  on  this  requirement,  and  replace  (4.7)  by  this 
condition. 

First,  we  compute  the  membership  function  of  the  fitted  function  F-.  Let 
Pt  be  the  panderance  matrix  of  the  fuzzy  model  parameter  vector  t.  (One 
obtains  an  estimate  of  Pt  concurrently  with  the  numerical  solution  of  the 
minimization  problem  (4.6),  see  Ref  3).  The  panderance  PFi  of  Fi  is  approxi¬ 
mately  given  by  Eq.  (3.27),  that  is,  by 

PFi  *  Fti  pt  F,i  •  <4-8) 

where  Fh  =  dF-/dt .  We  assume  that  the  r,  components  of  Ft  are  independent 
funtious  of  t  and  that  r,<p,  so  that  the  rank  of  F equals  r,.  Then  PFt  is  posi¬ 
tive  definite,  and  defines  the  following  norm  for  F,  : 

I  Fji  X  )||  pi  =  <piTpFilFi)I/2.  <4-9> 

and  the  corresponding  membership  function  in  the  A'-space 

FFj (X)  =  l-min|l,||Fi(X)||Fjj  .  (4.10) 

Next,  we  define  the  grade  of  compatibility  between  the  observation  A,  and  the 
fitted  model  by 


min jM|.(X),MFj(X)j  . 


(4. 11) 


The  constraint  (4.7)  of  the  minimization  problem  (4.6)  should  be  replaced  by 
the  requirement  that  the  grade  of  compatibility  is  for  each  observation  higher 
than  ,  that  is,  by 

XYC.  >  Y  ,  i  =  l,...,s,  (4.12) 

XFi 
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In  th  is  form  the  constraints  may  not  be  easily  implemented,  but  in  most 
cases  one  can  approximate  them  by  conditions  on  the  grade  of  collocation,  Eq. 
(3.31),  between  the  observations  A,  and  the  corrected  observations  A^+r,-. 
Using  Eq.  (3.32)  we  shall  express  these  conditions  in  terms  of  the  discord 
between  these  two  vectors.  We  use  ||c,  ||tt  as  a  distance  measure  h[Xi  +  ei  A^ 
from  A,  and  a  linearized  form  of  |  |E,  |  j/y  as  a  distance  measure  h(  A'.A^  +  r,) 
from  A, + c, .  Now  we  shall  derive  this  linearized  form  Let  Y0  be  a  point  on 
the  surface  F,  =  0.  Then  to  the  first  order 

AF|*  Fi(X0+b,t)-F|(Xo,»)*F,{X0+b.»)«FxI  b  (4.13) 

where  F^—dF^  A0  l)/dX.  Substituting  this  expression  for  E,  into  Eq.  (4.9) 
one  obtains 


Let 


PM  FrJ  PF?  Fxi 


(4. 14) 

(4.15) 


l\  1  is  a  n4  Xn, -matrix  and,  because  in  general  n, >r,  (unless  the  model  is  a 
point)  the  matrix  is  only  semidefinite  Therefore,  the  surfaces  b  TPh  lb  — const 
are  degenerated  ellipsoids  in  the  A’-spa.  e,  namely,  cylinders  with  rulings  that 
are  parallel  to  E,=  0.  The  pseudonorm  |  |H  If,  =  (  b  T  Fh 1  b)  only  measures 
that  component  of  b  which  is  orthogonal  to  Ft  =0.  (Orthogonal  in  the  sense  of 
the  inner  product  associated  with  the  pseudonorm).  The  use  of  the  linearized 
pseudonorm  IIMIf.  instead  of  ||AE,||f,  in  the  definition  (4.10)  of  the 
membership  function  ft f,  amounts  to  a  linearization  of  ft F ,  in  the  A-space. 
Now  we  use  IIMIf.  35  a  measure  for  the  distance  in  A'-space  from  the 
corrected  observation,  that  is. 


h(Xj  +  c,+b,  Xj+cj)  s  b  p •  .  (4.16) 


With  these  definitions,  the  discord  between  and  A, + c-  is  according  to  Eq 
(3.29)  r  . 


D;  5  min  max 


{NIxidMIn} 


(4.17) 
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where  c  =  A'-  A,  and  b  =X-  (A^+c,) 

If  the  linearization  (4.14)  is  sufficiently  accurate  and  the  dimension  of  Ft 
equals  one  then  the  value  of  the  discord  Z),  can  be  explicitly  calculated  from  the 
two  norms  of  the  residual  1 1 r,  j ) ,\7  and  ||c,  j|fi,  because  the  minimum  in  Eq. 
( 4.17)  is  attained  at  a  point  on  the  straight  line  segment  between  Xj  and  A,  and 
c,  .  We  demonstate  this  property  of  the  minimum  by  observing  that  c\~,  is  a  con¬ 
vex  function  and  showing  that  the  level  planes  =  const  are  tangential  to  the 
level  surfaces  of  c\;=  const  for  all  points  of  the  segment.  A  plane  by,  =  const  is 
spanned  by  vectors  q  which  satisfy  the  equation 

Fxi  Q  s  0  (4.18) 

(See  the  definition  F'q  ( 4.14).)  A  plane  through  the  point  A"  4-  c  and  tangen¬ 
tial  to  ||c  | }y;  =  const  is  spanned  by  vectors  q  which  satisfy  the  equation 

cVq  =  0  (4  .19  ) 

(See  the  definition  Eq  (4.5)).  Let  A'+r  be  on  the  straight  line  through  A,  and 
A’+c,  i.e  ,c  =  oc(.  Thn  js  a  solution  of  the  least  squares  problem  ( 4.6) 

Therefore,  it  satisfies  the  following  equation  from  Reference  4. 

ci  =  Pj^xi  (^xi^i^xi^)  ^xicj  •  (4.20) 

Substituting  o  times  the  right  hand  side  of  Eq.  (4.20)  for  r  in  Eq.  (4.19),  one 
obtains 


CiT  F*iT(FxiPiPKiTr‘  5  *  0  (4.21) 

The  set  of  all  vectors  q  satisfying  Eq  (4.21)  are  in  the  plane  tangential  to 
|  |r  j  |yj  =  const.  Because  q  satisfies  Eq  (  4.18),  it  also  satisfies  Eq.  (  4.21 )  and,  is 
therefore,  a  tangent  vector  Hence  1 16  |  |f,=  const  is  tangential  to 

||f  |  |y-=  const . 

Now  we  compute  the  discord  Pt.  taking  in  Eq  (4  17)  the  maximum  onl> 
over  points  of  the  straigt  line  between  A,  and  A^+c,.  On  that  line,  r  =  or,  and 
b={a  1 )  c,  Therefore,  Eq  (  4.17)  becomes 
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Oj  =  min  max 
a 

which  has  the  solution 


Di 


(4.22) 

(4.23) 


The  corresponding  grade  of  collocation  between  Aj-  and  Aj  +  c,  is 


X;  2  max 


{o,l-D,j  . 


(4.24) 


We  use  this  q,-  as  an  approximation  to  and  replace  the  compatibility 
condition  ( 4.12)  by 


Xj  >  r  ,  i  =  l  , . . . ,  s  . 


(4.25) 


With  the  same  approximation  we  define  the  overall  discord  between  data 
and  model  as 


D  =  mox  Dj  , 
i  =  l,  .,$ 


(4.26) 


and  the  overall  grade  of  compatibility  as 

y  =  rnin/j  =  max{o,l-D} 


(4.27) 


The  reformulated  least  squares  model  fitting  problem  (4  6)  and  (4.7)  is: 


rriinirrize  W  2  £  C  j  Tp/lq 

i«l 

subiect  to  ^(Xi  +  Ci  *  t)  =0.  i-1.2 


(4.28) 


and 


(4.29) 


The  numerical  solution  of  the  model  fitting  problem  may  be  obtained  by 
first  solving  Eq  (4  28).  and  checking  the  conditions  (4.29)  after  solution.  If 
th  e  conditions  are  satisfied,  then  one  has  the  solution.  If  for  some  A,  the  com¬ 
patibility  grade  is  less  than  the  required  q  .  then  one  may  take  several  actions, 
f  irst  one  should  recheek  the  data  and  depending  on  the  circumstances,  one 
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£ 

may  discard  outliers,  lower  the  required  compatibility  grade  7 ,  modify  the 
spread  estimates  of  observations  and/or  model,  change  the  model,  or  try  to 
find  a  solution  that  does  satisfy  Eq.  (4.29).  The  latter  means,  however,  an 
accommodation  of  outliers  at  the  expense  of  the  overall  fit.  A  solution  will  not 
always  exist  and,  if  it  does  exist,  it  usually  will  not  be  desirable.  Discarding  of 
observations  and  changing  the  model  equations  are  problem  dependent  actions 
and  will  not  be  discussed  here.  Changing  of  7  merely  changes  the  right  hand 
side  of  Eq.  (4.29).  The  more  general  actions  are  the  modifications  of  the 
spread  estimates  which  we  now  shall  discuss. 

We  first  consider  the  effect  of  changing  all  spread  estimates  of  the  obser¬ 
vations  by  a  factor  0  y  Then  all  panderance  matrices  Pt  are  multiplied  by 
0^-and  llr.ll/y  as  well  | |r,| |/r,  are  multiplied  by  1/0  A-.  (^1  in  Eq.  (4.8)  receives 
the  factor  0^-.  which  affects  ||6||f,  through  Eq.  (4.14)).  Consequently,  the 
discords  Di  in  Eq  (4.23)  are  multiplied  by  l/0y  too.  Hence  in  order  to  make 
all  discords  less  than  1  -  7  ,  i.e. ,  all  compatibility  grades  larger  than  7  ,  one  may 
multiply  the  estimated  data  spreads  by 

=  mQx|oi/{l  -  x)j.  (4.30) 

We  call  0 x  the  data  spread  dilator.  If  the  conditions  (4.29)  are  satisfied,  then 
the  data  spread  dilator  is  less  or  equal  one.  If  the  compatibility  grade  7 
between  data  and  model  is  less  than  7  .  then  the  data  spread  dilator  0 x  pro¬ 
vides  a  measure  for  the  seriousness  of  maximum  deviation  from  7  .  The  indi¬ 
vidual  0  vr  i.e.,  the  right  hand  side  terms  of  Eq.  (4.30),  provide  such  a  measure 
for  each  observation  A", 

Now  we  consider  model  spread  dilators.  One  may  use  such  a  dilator  if  the 
compatibility  grade  7  between  data  and  model  is  less  than  the  desired  7  .  but 
the  data  spread  estimates  and  the  model  equations  are  firmly  established.  Then 
one  may  achieve  the  desired  grade  of  compatibility  by  increasing  the  computed 
fuzziness  of  the  model,  i.e.,  by  dilating  the  spread  of  the  model  In  general, 
this  can  be  done  by  modeling  the  spreads  of  the  model  function  concurrently 
with  the  modeling  of  the  event.  We  consider  here  the  simplest  case  of  such  a 
modeling,  namely  the  multiplication  of  the  panderance  matrix  P,  in  Eq.  (4.8) 
by  a  factor  A  value  of  the  factor  may  be  obtained  from  one  of  the  following 
two  requirements.  First,  one  obtains  a  minimal  model  spread  dilator  0  p  by  the 
requirement  that  the  constraints  ( 4.29)  be  satisfied.  This  dilator  is  a  comple¬ 
ment  to  the  data  spread  dilator  0V.  We  obtain  its  value  by  observing  that  the 
multiplication  of  the  model  parameter  spreads  by  <?F  changes  ||c,||f-,  in  Eq 
(4.29)  by  the  factor  \/Cy.  and  leaves  the  ||c,||A7  unchanged.  Then  all 
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constraints  (4.29)  can  be  satisfied  if  <p F  has  the  value. 

*f'  m<ixH,ox|o,Ici|F./(i-;)-|ci||Fj/|ei|x,J  .  (4.31) 

qf  is  less  or  equal  one  if  all  £>,<1-'),  i.e.,  if  the  original  model  is  compatible 
with  the  data  to  the  desired  grade.  Otherwise  <j>F  has  a  value  bigger  than  one. 
The  individual  <pFi,  i.e.,  the  terms  on  the  right  hand  side  of  Eq.  (4.31)  again 
provide  a  measure  of  the  seriousness  of  the  incompatibility  between  the  obser¬ 
vation  and  the  fitted  model.  If  the  support  of  the  function  intersects 

with  the  surface  Ft  =  0  then  6 Fi  equals  zero. 

Another  model  spread  dilator  is  more  practical  for  application.  We  call  it 
the  inelusiie  model  spread  dilator  4> F  and  determine  its  value  from  the  require¬ 
ment  that  the  support  of  /i  7  should  include  all  support  ellipsoids  of  the 
functions  /jFi~*.  We  obtain  an  approximate  value  of  4> F  by  the  following 
algorithm.  First,  we  find  for  each  observation  A^,  that  point  on  the  straight  line 
through  A,  and  A,  +  c,  which  has  the  membership  value  7  and  which  is  on  the 
side  of  A,  opposite  to  A,  +  c,.  That  point  has  the  coordinates 

Yi  •  Xr«i(l-;)/|c,  ||xi  .  (4.32) 


The  membership  value  of  /,  at  that  point  is 

^F.(Yi)  *  1-min jlf||Fj(Yj)|FlJ  ,  (4.33) 

where 


II  F(  (v i )  [F j  -  (fs (r j )TPFri  f,  (v, )) 1/2 


is  calculated  by  Eq  (4  9).  We  want  this  value  to  be  equal  or  larger  than  7  .  Therefore, 
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the  spread  of  the  model  parameter  vector  should  be  multiplie  d  by 


BFi(vi )||Fj / O-r)  • 


(4.34) 


The  overall  inclusive  model  spread  dilator  is 

<X>F  =  max  <*>Fj  .  (4.35) 

i  ri 

In  summary,  we  have  introduced  the  following  indicators  which  measure  the 
quality  of  the  model  fitting: 

D  -  d  iscord  between  data  and  model, 
q  -  grade  of  compatibility  between  data  and  model, 

0  v  -  data  spread  dilator, 

<>  f  -  minimal  model  spread  dilator, 

<t 'p  ~  inclusive  model  spread  dilator. 

The  last  three  indicators  depend  on  a  desired  minimum  grade  of  compatibility, 
whereas  the  first  two  indicators  are  independent  of  such  a  grade. 


5.  EXAMPLE 


As  an  example  of  a  model  with  fuzzy  experimental  data  we  present  the  fol¬ 
lowing  treatment  of  perforation  of  a  plate  by  projectiles.  In  perforation  experi¬ 
ments,  one  is  interested  in  the  minimum  striking  velocity  which  is  needed  to 
achieve  perforation.  That  velocity  is  called  the  ballistic  limit  velocity  and  for  a 
given  target  it  depends  on  the  materials  of  the  projectile  and  target,  on  the 
form  of  the  projectile,  and  on  the  obliquity  of  the  impact.  The  limit  velocity 
typically  is  determined  by  firing  projectiles  with  different  velocities  and 
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observing  their  residua]  velocities  (the  velocity  alter  penetration).  The  limit 
velocity  is  computed  from  these  data  by  one  of  several  available  methods  which 
all  amount  to  some  kind  of  interpolations  between  zero  and  non-zero  residual 
velocities.  The  result  contains  observational  errors  as  well  as  intrinsic  varia¬ 
tions  of  the  limit  velocity  between  rounds  due  to  unknown  causes,  e.g.,  due  to 
variations  of  target  and  projectile  material  properties.  (Such  variations  manifest 
themselves  e.g.,  by  overlapping  regions  of  zero  and  non-zero  residual  veloci¬ 
ties.)  The  number  of  experiments  usually  is  too  small  to  determine  a  distribu¬ 
tion  of  the  result,  but  expert  estimates  of  possible  variations  are  available  from 
past  experience.  Table  I  contains  a  set  of  observed  limit  velocities.  (The  data 
are  taken  from  Ref.  5).  The  targets  were  identical  in  all  cases.  The  projectiles 
were  cylindrical  rods  with  different  nose  shapes,  quantitated  in  the  table  by  a 
"shape  quotient”  which  is  the  ratio  of  nose  length  to  rod  length.  Because  it  is 
not  certain  that  this  ratio  adequately  describes  the  significance  of  the  projectile’s 
geometry,  a  substantial  spread  of  0.05  was  assigned  to  the  value  of  the  ratio. 
The  spread  of  the  limit  velocity  was  estimated  from  past  experience  to  be  about 
50  m/s.  However,  in  two  cases  the  actual  observations  were  not  sufficient  to 
guarantee  such  an  accuracy,  and,  therefore,  corresponding  larger  spreads  were 
assumed  in  these  cases. 

The  projectiles  were  fired  in  two  modes:  normal  to  the  target  plate,  and  at 
a  constant  angle  of  obliquity.  The  purpose  of  the  experiments  was  to  deter¬ 
mine  whether  and  how  the  shape  of  the  projectile  affects  the  limit  velocity. 
The  problem  was  handled  by  fitting  a  constant  and  a  linear  function  to  the  two 
series  of  observed  limit  velocities,  and  by  comparing  the  results.  The  numeri¬ 
cal  fitting  of  the  fuzzy  data  was  done  with  the  general  least  squares  computer 
program  COLSAC 

The  data  are  in  the  present  case  two-component  fuzzy  vectors  ( Q ,  V). 
The  component  Q  is  the  dimensionless  shape  quotient  and  the  component  V  is 
the  corresponding  limit  velocity  in  m/s.  The  models  of  the  limit  velocity  were 
the  linear  functions 


V  =  a  +  bQ  . 


(5.1) 


L’sing  the  technique  described  in  Section  4,  one  obtains  as  a  result  of  each 
model  fitting  a  fuzzy  parameter  vector  with  the  components  a  (m/s)  and  b 
(m/s),  and  a  corresponding  (2X2)  panderance  matrix.  In  case  of  a  constant 
model  one  sets  6=0  in  Eq.  (5.1)  (6  is  a  crisp  zero  in  that  case),  and  the  model 
fitting  algorithm  provides  only  a  and  its  panderance.  In  the  (a,6)-space,  the 
special  solution  vector  (a,0)  with  a  crisp  second  component  has  a  singular 
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panderance  matrix.  The  support  of  the  corresponding  membership  function  is  a 
segment  of  the  a-axis.  The  numerical  values  of  the  four  fuzzy  parameter  vec¬ 
tors  are  listed  in  Table  II.  The  data  and  the  fitted  curves  are  shown  in  Figures  2 
through  7. 

Next,  we  discuss  and  illustrate  the  various  measures  of  compatibility 
between  data  and  models.  Table  II  lists  these  indicators  for  each  data  set 
(  Q,  V)  from  the  normal  impact  observations  series.  (See  also  Figures  2  and  4). 
To  be  specific  we  have  assumed  that  the  desired  grade  of  compatibility  is 
7  =0.3,  and  Table  III  shows  that  the  normal  impact  data  exceed  this  level  of 
compatibility  with  the  constant  model  as  well  as  with  the  linear  model.  Conse¬ 
quently  the  data  dilators  and  the  minimal  model  dilators  all  are  less  than  one. 

Figure  2  shows  the  normal  impact  data  and  the  fitted  constant  model. 
Each  data  point  is  plotted  with  the  support  ellipse  of  its  membership  function. 
A  dotted  ellipse  indicates  the  level  line  pi  x  =  0.3.  The  fitted  velocity  model  is 
the  line  Y=  810.2  m/s.  It  is  a  fuzzy  line,  with  a  spread  of  24.5  m/s,  as  indi¬ 
cated  in  Figure  2.  The  dotted  horizontal  lines  are  the  level  lines  pi x  =  0.3. 
The  data  dilators  and  the  minimal  model  dilators  all  are  less  than  one  (see 
Table  III)  because  the  supports  of  pi  ^ 7- 0.3  and  pi  F~  0.3  intersect.  The  inclusive 
model  dilators  in  Table  III  for  the  desired  grade  of  compatibility  7  =  0.3  are 
larger  than  one  because  the  supports  of  pi  a?- 0.3  are  not  subsets  of  the  support 
pi  f-  0.3.  Similarly,  the  supports  of  pix-  are  not  subsets  of  the  support  of  piF, 
and  therefore,  also  the  inclusive  model  dilators  for  7  =0  ought  to  be  larger 
than  one.  Their  maximum  value  <F0  =  6.832  is  listed  in  Table  II.  If  the  spread 
of  the  model  is  dilated  by  the  factors  <F0  or  $03,  then  it  includes  the  supports 
of  the  data  up  to  the  corresponding  level  0.0  or  0.3,  respectively. 

Figure  3  shows  the  same  data  as  Figure  2  but  the  fitted  model  is  a  linear 
function.  The  discords  between  data  and  fitted  model  are  smaller  than  for  the 
constant  model  (see  Table  III).  The  inclusive  model  dilators  have  reasonable 
values,  4>0  =  2.290  for  7  =  0.0  and  4>0  3  =  2.407  for  7  =  0.3.  Figure  4 
shows  the  effect  of  these  dilators.  It  displays  the  same  data  and  model  as  Fig¬ 
ure  3,  but  the  model  spread  has  been  dilated.  The  solid  line  indicates  the 
boundary  of  the  support  of  pi  F,  if  the  dilator  4>0  =  2.290  is  used.  The  dotted 
lines  indicate  the  boundary  of  the  support  of  piF  and  the  level  line  pi  F  =  0.3  if 
the  dilator  <t>0  3=2.407  is  used.  In  the  former  case,  the  support  of  pi F  includes 
all  supports  of  the  data.  In  the  second  case,  the  support  of  pi  F-  0.3  includes  the 
supports  of  pi xr  0.3.  Figure  4  confirms  that  the  approximately  calculated  dila¬ 
tors  by  Eq.  (4.34)  are  sufficiently  accurate  in  this  example. 

Next,  we  consider  the  oblique  impact  data  listed  in  Table  FV.  The  table 
indicates  that  the  data  are  incompatible  with  a  constant  limit  velocity  model, 
and  Figure  5  shows  that  the  supports  of  the  data  indeed  do  not  intersect  with 
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the  support  of  the  fitted  model,  except  for  one  data  point.  All  dilators  are 
larger  than  one,  particularly  the  inclusive  model  dilators  $0=11.579  and 
$0  3=14.941.  (Table  II). 

Figure  6  shows  the  result  of  fitting  a  linear  function  to  the  oblique  impact 
data.  The  data  set  is  compatible  with  the  model,  the  grade  of  compatibility  is 
0.298.  Consequently,  all  data  dilators  and  minimal  model  dilators  for  the 
desired  grade  of  compatibility  q=0.3  are  less  than  one  or  close  to  one  (Table 
IV).  The  inclusive  model  dilators  are  somewhat  large,  $0  =  4.625  and 
$0  3=5.464.  Using  these  dilators  to  increase  the  spread  of  the  fitted  function 
one  obtains  the  result  shown  in  Figure  7.  As  in  Figure  4,  we  have  plotted  as 
solid  curves  the  level  lines  Up  =  0,  dilated  by  $0'  F°r  the  spread  of  the  fitted 
function,  dilated  by  $03  we  have  plotted  as  dotted  curves  the  level  lines  up  = 
0.0  and  // p  =  0.3.  The  regions  between  these  curves  should  include  all 
corresponding  supports  of  the  data  membership  functions.  One  notices  a  slight 
intersection  of  the  boundary  curves  of  the  second  data  point  ( <2=0.052)  with 
the  model  spread  boundaries.  This  intersection  (instead  of  a  tangent)  is  caused 
by  the  approximate  calculation  of  the  dilators  by  Eq.  (4.34),  only  using  function 
values  along  the  line  through  the  observation  A'  and  the  corrected  observation 
A'+c. 

One  can  conclude  from  the  presented  results  that  in  case  of  oblique  impact 
the  limit  velocity  indeed  depends  on  the  nose  shape,  increasing  with  the 
slenderness  of  the  projectile's  nose.  In  the  normal  impact  case  there  seems  to 
be  a  slight  decrease  of  the  limit  velocity  with  increasing  nose  shape  quotient. 
The  data  are  compatible  with  the  constant  model  to  a  grade  0.640  and  compati¬ 
ble  with  the  linear  model  and  a  grade  0.859.  The  increase  of  the  compatibility 
is  not  very  large  and  one  might  ask  to  what  degree  both  fitted  models  differ 
from  each  other.  A  measure  for  this  difference  is  the  grade  of  collocation 
between  the  two  fuzzy  model  parameter  vectors  which  define  the  models 
(Table  II).  Figure  8  shows  the  two  vectors  with  the  corresponding  supports  of 
their  membership  functions  in  the  parameter  space.  As  mentioned  before,  the 
support  of  the  parameter  (  a ,0) ,  i.e . ,  the  constant  model  parameter,  consists  of 
a  segment  of  the  a-axis.  The  support  of  the  linear  model  parameter  member¬ 
ship  function  is  aD  ellipse,  shown  in  Figure  8  with  a  solid  line.  The  dotted  line 
curve  in  the  figure  is  the  level  line  /i  =  0.3.  The  grade  of  collocation  between 
these  two  vectors  turns  out  to  be  0.080.  If  one  increases  the  spread  of  the 
parameters  by  the  factors  $0  or  <t>0  3,  then  one  also  increases  the  grade  of  collo¬ 
cation,  namely  to  0.598  and  0.618,  respectively. 

For  the  oblique  impact  data  the  model  vectors  of  constant  and  linear 
models  have  a  zero  grade  of  collocation.  Figure  9  shows  the  situation  of  the 
vectors  in  the  parameter  space.  The  discord  between  the  vectors  is  in  this  case 
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3.250.  (If  one  uses  the  dilators  4>0  or  4>03,  then  the  grades  of  collocation  are 
0.296  and  0.404,  respectively.)  The  large  discord  is  an  indication  that  there  is  a 
significant  difference  between  both  models.  In  such  a  situation  one  should  use 
the  better  model  to  represent  the  data,  in  this  case  the  linear  model.  This  leads 
to  the  same  conclusion  as  above:  for  oblique  impact,  the  limit  velocity  does 
increase  with  the  slenderness  of  the  projectile. 

We  make  a  final  comment  to  the  representation  of  the  data  by  a  model, 
and  consider  as  an  example  the  oblique  impact  data.  A  straight  forward  deter¬ 
mination  of  the  model  produces  the  result,  shown  in  Figure  6,  which  is  compa¬ 
tible  with  the  data,  but  does  not  include  all  supports  of  the  data.  Using  model 
dilators  one  obtains  the  results  shown  in  Figure  7,  which  includes  all  data  sup¬ 
ports,  but  seems  to  be  excessively  pessimistic  at  some  parts  of  the  fitted  curve. 
Hence,  if  one  intends  to  use  the  fitted  model  for  predictive  purposes,  one  may 
also  want  to  model  the  spread  of  the  results  so  that  it  reflects  the  available  data 
more  closely.  Such  a  modeling  depends  very  much  on  the  particular  applica¬ 
tion.  In  the  present  example,  using  a  spread  model  that  is  of  second  degree  in 
Q ,  and  symmetric  about  the  fitted  line  one  obtains  the  result  shown  in  Figure 
10.  It  is  probably  more  reasonable  for  predictive  purposes  than  the  result 
shown  in  Figure  7. 


6.  SUMMARY  AND  CONCLUSIONS 


This  article  exposes  several  possible  fuzzifications  of  model  fitting,  and 
treats  in  detail  the  fitting  of  a  crisp  model  to  fuzzy  vector  data.  The  result  is  a 
fuzzy  fitted  model,  whereby  the  fuzziness  of  the  model  depends  of  the  fuzzi¬ 
ness  of  the  data. 

The  data  type  is  specialized  to  fuzzy  vectors  with  conical  membership  func¬ 
tions.  The  restriction  to  conical  functions  greatly  facilitates  the  numerical  treat¬ 
ment  of  the  problem,  and  a  simple  spread  propagation  formula  can  be  derived 
for  such  functions.  A  new  concept  is  the  discord  between  two  fuzzy  vectors.  It 
is  a  dimensionless  and  symmetric  measure  of  the  separation  of  the  vectors.  A 
complementary  new  concept  is  the  grade  of  collocation  of  two  elements  of  a 
fuzzy  vector  space. 

The  restriction  to  conical  membership  functions  permits  one  to  use 
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available  model  fitting  computer  programs  with  little  changes.  The  quality-  of 
the  fit  is  characterized  by  a  discord  between  data  and  fitted  model  and  a  grade 
of  compatibility  between  data  and  model.  Both  concepts  naturally  arise  from 
the  problem  formulation.  If  the  grade  of  compatibility  is  not  as  high  as  desired, 
then  one  can  either  try  a  different  model,  or  increase  the  fuzziness  of  the  data 
and/or  of  the  fitted  model,  such  that  it  becomes  compatible  with  the  data.  The 
factors  which  achieve  this  are  defined  as  data  and  model  spread  dilators,  respec¬ 
tively,  and  explicit  approximate  formulas  are  derived  for  their  calculation.  The 
values  of  the  dilators  depend  on  the  desired  grade  of  compatibility,  as  well  as 
on  the  data  and  their  membership  functions. 

The  application  of  model  fitting  to  fuzzy  data  is  likely  to  be  in  problems 
with  scarce  data  and  little  information  about  data  distribution.  The  description 
of  the  observed  event  by  a  fuzzy  function,  as  derived  in  this  paper,  explicitly 
exposes  the  lack  of  exact  information.  A  proper  quantitation  of  this  deficiency 
is  achieved  by  the  use  of  model  dilators.  The  restriction  of  the  algorithms  to 
conical  membership  functions  is  likely  not  a  serious  impediment:  because  of 
the  simplicity  of  the  algorithms  one  easily  can  carry  out  the  model  fitting  with 
different  estimates  of  data  spreads  and,  comparing  the  results,  determine  the 
significance  of  such  estimates. 

In  summary,  the  described  model  fitting  provides  realistic  representation  of 
events  about  which  only  fuzzy  information  is  available. 
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Projecti  le 
Type 


Flat 

Hemi spheri  ca 
Snort  Cone 
Long  Cone 
Full  Cone 


Table  I.  Ballistic 


Normal 


Shape 

Limi  t 

Quotient 

Spread 

Velocity 

Q 

SQ 

V  (m/s) 

0.000 

0.050 

841 

0.052 

0.050 

846 

0.165 

0.050 

803 

0.232 

0.050 

795 

0.431 

0.050 

772 

Limit  Data 


Impact 

Oblique 

Limit 

Impact 

Spread 

Velocity 

Spread 

sv  (m/s) 

V  (m/s) 

sv  (m/s) 

50 

1060 

50 

50 

978 

90 

50 

1111 

50 

50 

1255 

50 

120 

1300 

50 

t 


i 


I 

i 

1 


2" 


A 


Table  II.  Parameters  of  Fitted  Models 


Model  V  =  a  +  bQ 
Normal  Impact 


Model  Type 

a 

s 

a 

b 

s 

b 

c 

ab 

<t> 

0 

<t> 

0.3 

Constant 

819.2 

24.5 

0.0 

0.0 

0.0 

6.832 

7.658 

Linear 

845.3 

37.9 

-207.8 

226.0 

-0.752526 

2.290 

2.407 

Oblique  Impact 

Model  Type  a 

sa 

b 

Sb 

c 

ab 

*0 

0.3 

Constant 

1166.9 

24.1 

0.0 

0.0 

0.0 

11.579 

14.941 

Li  near 

1030.2 

51.2 

688.7 

211.4 

-0.824466 

4.625 

5.464 

The  dilators  $  are  not  included  in  the  quoted  values  of  the  spreads  s 
and  sb.  The  c  ^  are  the  concordances  between  the  components  a  and  b  of 
the  parameter  Actors.  The  units  of  a,  s,,  b  and  sK  are  m/s. 
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Table  111.  Model  Fitting  Results  for  Normal  Impact  Data  Constant  Model 


Constant  Model 


Projecti le 

| 

'c:* 

Discord 

Flat 

0. 

,43b 

0 

.891 

0. 

293 

Hemi spneri cal 

0. 

.536 

1 

.095 

0. 

360 

Short  Cone 

0. 

,324 

0 

.662 

0. 

218 

Long  Cone 

0. 

.484 

0 

.989 

0. 

325 

Full  Cone 

0. 

,393 

1 

.928 

0. 

327 

Linear  Model 


Projecti 1 e 

c 

"X 

C"F 

Discord 

Flat 

0. 

,084 

0. 

.114 

0. 

,048 

Hemi spheri cal 

0. 

,225 

0. 

.378 

0. 

,141 

Short  Cone 

0. 

157 

0. 

.302 

0. 

,103 

Long  Cone 

0. 

041 

0. 

,061 

0. 

,025 

Full  Cone 

0. 

13b 

U. 

,222 

0. 

.084 

Tie  dilators  are  calculated  for  a 


Comp. 

Data 

Model 

Di 1 ators 

grade 

Di  1 ator 

Minimal 

Incl usi ve 

0.707 

0.418 

0.0 

3.315 

0.640 

0.514 

0.0 

3.607 

0.789 

0.311 

0.0 

2.989 

0.675 

0.464 

0.0 

3.456 

0.673 

0.467 

0.0 

7.658 

Comp. 

Data 

Model 

Di  1  ators 

grade 

Di  1 ator 

Mi nimal 

Incl usi ve 

0.952 

0.069 

0.0 

1.461 

0.859 

0.202 

0.0 

2.309 

0.897 

0.148 

0.0 

2.407 

0.975 

0.035 

0.0 

1.615 

0.916 

0.120 

0.0 

1.944 

desired  compatibility  grade  of  0.3. 
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Table  IV.  Model  Fitting  Results  for  Oblique  Impact  Data 


Constant  Model 


Project i 1 e 

i  c . , 

A 

||C  ip 

Discord 

Flat 

2.138 

4.437 

1.443 

Hemi spher i ca 1 

2.099 

7.843 

1.656 

Short  Cone 

1.118 

2.322 

0.755 

cony  Cone 

1.762 

3.656 

1.189 

Full  Cone 

2.662 

5.525 

1.796 

Linear  Model 


Pro. 

ject  i  1  e 

’i 

,r  V 

Dis 

cord 

Flat 

0, 

.  4  u  1 

0. 

.611 

0. 

272 

Hem  s 

phen  cal 

0. 

.913 

1, 

.950 

0. 

622 

Short 

Cone 

0. 

.541 

1 . 

.066 

0. 

359 

Lone 

Cone 

1. 

.071 

2. 

.041 

0. 

702 

Full 

Cone 

0 . 

.445 

0. 

.495 

O. 

235 

Comp. 

Data 

Model 

Di 1 ators 

grade 

Di  1 ator 

Minimal 

Incl usi ve 

0.0 

2.062 

4.265 

8.417 

0.0 

2.366 

7.468 

14.94] 

0.245 

1.078 

1.241 

5.392 

0.0 

1.698 

3.148 

7.299 

0.0 

2.566 

5.817 

9.968 

Comp . 

Data 

Model 

Di 1 ators 

grade 

Di  1 ator 

Minimal 

Incl usi ve 

0.728 

0.380 

0.0 

1.889 

0.378 

0.905 

0.690 

5.464 

0.641 

0.518 

0.0 

3.697 

0.298 

1.027 

1.081 

5.283 

0.765 

0.328 

0.0 

1.641 

Tne  dilators  are  calculated  for  a  desired  compatibility  grade  of  0.3. 
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Inpact  Data 


VELOCITY .  M/S 


LIMIT  VELOCITY.  M/S 


SHAPE  QUOTIENT 


Figure  7.  Linear  Model  with  Dilated  Spread  Fitted  to  Oblique  Impact  Dat 
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Figure  9.  Model  Parameters  for  Oblioue  Impact  Data 
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